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摘要 :根据 毛乌素 沙 地 典型 站 点 近 60 年 (1955 一 2014 年 ) 逐日 气象 数据 ,利用 Penman-Monteith 公式 计算 毛乌素 沙 地 各 气象 站 
点 参考 作物 蒸 散 量 (reference crop evapotranspiration，E7 ) 及 研究 区 域内 整体 ETu。Mann-Kendall 趋势 检验 法 和 ArcGIS 的 协同 
克 里 格 插值 法 用 于 分 析 E7, 时 空 变化 特征 ,同时 ,利用 敏感 性 分 析 方 法 对 ET 的 变化 成 因 进行 分 析 。 结 果 表 明 :(1) 近 60 年 毛 
乌 素 沙 地 E7, 多 年 平均 值 为 1048.81mm ,年 际 变 化 呈现 缓慢 上 升 趋势 。 年 内 变化 夏季 最 高 ,冬季 最 低 。 区 域内 ET 空间 分 布 整 
体 呈 现 自 西 向 东 递 减 趋 势 。(2) 7, 年 变化 对 风速 的 敏感 程度 最 大 ,日 照 时 数 和 气温 次 之 ,相对 湿度 最 小 。 春 、 秋 两 季 PT 变化 
对 日 照 时 数 最 为 敏感 ; 夏 冬 两 季 E7, 变 化 对 相对 湿度 最 为 敏感 。 空 间 分 布 上 ,毛乌素 沙 地 东南 部 地 区 为 气温 敏感 系数 高 值 区 
西北 部 地 区 为 相对 湿度 和 日 照 时 数 敏 感 系数 高 值 区 ,南部 为 风速 敏感 系数 高 值 区 。(3) 通 过 计算 气象 因子 对 E7, 变 化 的 贡献 量 
得 出 ,气温 是 影响 毛乌素 沙 地 E7, 年 变化 的 主导 因子 。 夏 季 E7, 变 化 的 主导 因子 是 风速 ;春秋 、 冬 三 季 主 导 因 子 是 气温 。 空 间 
分 布 上 ,毛乌素 沙 地 西南 部 地 区 E7, 变 化 的 主导 因子 为 风速 ,东部 地 区 主导 因子 为 气温 。 

关键 词 :协同 克 里 格 插值 法 ;毛乌素 沙 地 ;敏感 性 分 析 ;Penman-Monteith 公式 ; Mann-Kendall 趋势 检验 法 
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Abstract: Reference crop evapotranspiration ( ET,) is often used to evaluate the capacity of atmospheric evaporation in the 
terrestrial hydrology cycle. Simulations with climate models and atmospheric general circulation models have shown that 
changes in meteorological factors affect £7T,. The Mu Us desert is located in the arid and semi-arid region of northwest 
China, which experiences severe water shortages. To better understand the ET, spatio-temporal distribution characteristics in 
the Mu Us desert, the Penman-Monteith equation recommended by the FAO in 1998 was used to calculate ET, at typical 
meteorological stations, as well as the overall Mu Us desert ET,, based on 60 years of daily meteorological data (1955— 
2014). A Mann-Kendall trend test was used to determine the ET spatio-temporal distribution characteristics. Furthermore, 
a Co-Kriging interpolation was used to draw ET, maps, which presented a clear spatial distribution of 有 Tu in the Mu Us 
desert. In addition, a sensitivity analysis was used to analyze the causes of kT, change. Annual average ET, was 1048.81 
mm for the past 60 years. There was a slight increasing trend in the inter-annual variation. The trend of annual £7, showed a 


single peak curve, with the highest value in June (5.77 mm/d), and the lowest value was in January (0.54 mm/d). The 
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average summer ET, was 5.23 mm/d, and average winter FET, was 0.70 mm/d. Spatially, ET, decreased gradually from west 
to east, ranging from 1058.53 to 1103.85 mm in the west, and from 1014.05 to 1043.44 mm in the east. The sensitivity 
analysis indicated that ET, was positively related to wind speed, air temperature, and solar duration, but negatively related 
to relative humidity. Furthermore, ET, was most sensitive to changes in wind speed, then solar duration, air temperature, 
and relative humidity, respectively. Seasonally, ET, was more sensitive to solar duration in spring and autumn, but more 
sensitive to relative humidity in summer and winter. Wind speed was the dominant controlling factor on ET, in summer, but 
in spring, autumn and winter, it was air temperature. Overall, air temperature was the dominant controlling factor on annual 
ET, in Mu Us desert. Spatially, ET, was more sensitive to temperature in the southeast region of the Mu Us desert, to 
relative humidity and solar duration in the northwest region, and to wind speed in the south. The spatial distribution of 
ET,'s dominant controlling factor was different from its sensitive coefficient. Wind speed was the dominant controlling factor 
on 上 ET，distribution in southwest of Mu Us desert, and air temperature in east. The meteorological factors here were 
considered dominant in controlling the spatio-temporal variations of the Mu Us desert 五 7,. However, advection variation 
caused by wind speed, temperature and relative humidity, inter-annual variation of the atmospheric circulation, surface 
uniformity, and human activities, also influence ET, in China’ s arid and semi-arid region. Therefore, the results presented 
here only reflect the meteorological variables that were included in the analyses. More comprehensive hydrological studies 


are therefore required in the future. 
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trend test 


参考 作物 蒸 散 量 (reference crop evapotranspiration, £7, ) 作为 水 文 循环 过 程 的 重要 参数 ,能 够 为 评价 气候 
干旱 程度 、 估 算 作物 需 水 量 、 生 产 潜力 以 及 水 资源 供需 平衡 等 提供 依据 ,是 农业 生产 中 实现 合理 调配 水 土 
资源 及 缓解 水 资源 亏 缺 的 重要 参考 值 。 气 候 模式 和 大 气 环流 模式 的 数值 试验 表明 ,气候 因素 的 变动 将 直接 影 
响 ET 的 变化 。 在 气候 变化 背景 下 ,气温 、 相 对 湿度 太阳 辐射 等 气象 因子 的 波动 都 将 导致 57, 不 同 程度 
的 变化 。 因 此 ,分析 ET 的 变化 特征 及 其 变化 成 因 , 有 助 于 深入 理解 区 域 水 循环 对 气候 变化 的 啊 应 机 制 。 

关于 £7 的 研究 在 全 球 范围 内 已 广泛 开展 ,北半球 太阳 辐射 减少 ,南半球 云 量 、 气 溶胶 浓度 增加 均 被 认为 
是 全 球 E7, 减 少 的 主要 原因 '**i 。 同 时 ,不 同 地 区 87, 变化 对 气象 因子 的 响应 程度 有 所 差异 ,湿润 半 湿 润 地 区 
ET 变化 与 太阳 辐射 关系 密切 ,而 风速 ,温度 和 相对 湿度 对 干旱 半 干 旱地 区 E7, 变 化 作用 更 为 显著 '”" 。 过 去 
数 十 年 中 国 7, 总 体 呈 下 降 趋 势 外 ,但 东北 地 区 松 嫩 平原 西部 和 黄河 上 游 BT 又 呈 微 弱 上 升 趋势 "5" , 西 
北 地 区 ET 人 研究 结果 差异 较 大 ,乌鲁木齐 河流 域 及 黄土 高 原 以 西 地 区 ET 呈现 下 降 趋势 ,但 是 西 辽河 流域 及 黄 
土 高 原 以 东 ET 呈现 上 升 趋势 '“” 。 由 于 E87, 的 变化 受到 气象 要 素 、 地 形 .植被 .研究 区 域 和 时 段 等 多 因素 共 
同 影响 ,所 以 局 地 的 ET 变化 可 能 不 同 ,不 同 研 究 结果 存在 偏差 。 现 有 研究 表明 ,ET 的 变化 特征 及 其 主导 
子 具 有 明显 的 区 域 性 差异 。 基 于 1957 一 2012 年 全 国 气象 资料 分 析 得 出 气温 是 影响 中 国 57, 的 主要 因子 ”|， 
而 西北 地 区 ,华北 平原 和 川中 丘陵 B7, 变 化 的 主要 归 因 于 风速 和 日 照 时 数 '** 。 

毛乌素 沙 地 是 中 国 四 大 沙 地 之 一 ,位 于 中 国 西北 部 ,行政 区 地 跨 内 蒙古 自治 区 、 陕 西 省 和 宁夏 回族 自治 
区 ,降水 波动 性 大 ,旱季 年 份 居 多 ,水 资源 短缺 是 限制 毛乌素 沙 地 生态 环境 恢复 和 农业 生产 潜力 的 重要 因 
素 '*””] 。 随 着 全 球 气候 变化 程度 加 剧 ,毛乌素 沙 地 气候 呈现 暖 干 化 趋势 ,ET, 也 随 之 发 生变 化 :2 ,但 该 地 区 
ET 变化 的 主导 因子 尚 不 明确 。 本 文 基于 毛乌素 沙 地 典型 气象 站 点 1955 一 2014 年 气象 资料 ,利用 Penman- 
Monteith 公式 计算 毛乌素 沙 地 ET ,通过 Mann-Kendall 趋势 检验 法 和 ArcGis 的 协同 克 里 格 插值 法 ,探究 毛 乌 
素 沙 地 近 60 年 ET 时空 变化 特征 ,并且 利用 敏感 性 分 析 法 定量 分 析 各 气象 因子 对 E7, 的 影响 作用 。 研 究 结 有 果 
有 助 于 深入 理解 本 区 域内 气候 变化 对 水 循环 的 影响 ,同时 为 该 区 农业 生产 中 水 土 资源 的 合理 配置 提供 指导 。 
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1 研究 区 概况 


毛乌素 沙 地 位 于 中 国 西北 部 干旱 半 干 旱地 区 ,地 理 位 置 介 于 37*30 一 39°20'N 和 107°20' 一 111?30 中 之 
间 , 包括 内 蒙古 自治 区 鄂尔多斯 南部 ,陕西 省 榆林 地 区 北部 ,以 及 宁夏 回族 自治 区 盐池 县 东北 部 , 沙 地 面积 约 
4 万 km 。 本 区 域 是 鄂尔多斯 高 原 向 陕 北 高 原 过 渡 地 带 , 地 势 呈现 由 西北 向 东南 倾斜 ,属于 温带 季风 区 ,年 平 
均 气 温 为 6.10 一 8.15 % ,降水 量 呈 波动 性 变化 。 区 域内 自然 环境 恶劣 ,生态 和 人 口 的 双重 压力 严重 影响 当地 
农业 生产 和 社会 经 济 发 展 。 


2 研究 方法 


2.1 数据 来 源 

本 文选 取 毛 乌 素 沙 地 内 部 和 周边 8 个 气象 台 站 气象 资料 ,其 中 ,内 蒙古 自治 区 3 个 (鄂尔多斯 市 、 鄂 托 克 
旗 和 伊 金 截 治 旗 ) ,陕西 省 4 个 ( 定 边 县 横山 县 榆林 市 和 靖 边 县 ) , 末 夏 回族 回族 自治 区 1 个 (盐池 县 ) , 靖 边 
县 站 气象 资料 由 于 年 代 序 列 较 短 , 仅 用 于 分 析 毛 乌 素 沙 地 参考 作物 蒸 散 量 空间 分 布 特征 。 和 气象 资 料 来 源 于 国 
家 气象 信息 中 心 (http://www.escience.gov.cn/) ,为 1955 一 2014 年 逐日 气象 变量 ,包括 最 高 气温 .最低 气温 . 平 
均 气 温 、 相 对 湿度 .平均 风速 和 日 照 时 数 。 采 用 气象 学 标准 进行 季节 划分 , 即 春季 为 3 一 5 月 ,夏季 为 6 一 8 月 ， 
秋季 为 9 一 11 月 ,冬季 为 12 一 次 年 2 月 。 
2.2 ”Penman-Monteith 公式 计算 

本 文中 参考 作物 蒸 散 量 (ET ,也 称 潜在 蒸 散 量 ) 计算 采用 世界 粮农 组 织 (FAO ) 推荐 并 修订 的 Penman- 
Monteith 公式 。 表 达 式 为 : 


900 
0.408A(R, -G6G) +y F732. — e,) 


ET, = 1 
A+y(l + 0.34U,) (1) 


式 中 ,E7, 为 参考 作物 燕 散 量 (mm/d) ;RR, 为 参考 作物 冠 层 表面 净 辐 射 量 (MJ (md ) ;6 为 土壤 热 通 量 (MJ 
md !) ;7 为 平均 气温 (°C ) ;ie 为 饱和 水 汽 压 (kPa) ;e, 为 实际 水 汽 压 (kPa) ;A 为 水 汽 压 对 温度 的 斜率 (kPa/ 
% ) ;7 为 干 湿 球 常数 (kPa/5C ) ;U, 为 距 地 面 2m 高 处 风速 (m/s) 。 

Penman-Monteith 公式 中 辐射 项 采用 经 辐射 校正 的 FA056-PM 模型 计算 '” ,表达 式 为 ; 

R,=(1- [02 + 079 轴 ]R -| (056 = 0.25 /6-) [09 + oj (2) 
NS 2 . A N 

式 中 ,o 为 波 尔 兹 曼 常 数 [4.903x10”MJA(K* m d) ];R. 为 晴天 辐射 (MJ m”d ') ;n 为 实际 日 照 时 数 (h),N 
为 可 照 时 数 (h) ,a 为 地 表 反 射 度 , 取 值 0.23。 
2.3 气象 因子 敏感 性 分 析 

参考 作物 蒸 散 量 对 气象 因子 的 敏感 系数 计算 :2 , 即 通过 计算 FI7, 对 各 气象 因子 变量 的 偏 导数 ,将 气象 因 
子 对 E7, 的 变化 无 量 纲 化 ,定量 表达 BT, 对 单个 气象 因子 的 敏感 性 。 表 达 式 为 ， 
AET /ET,) 9ET, yi 
[ AD | ow ET, 3) 
式 中 ,5, 为 ET, 对 气象 因子 变量 的 敏感 系数 ,无 量 纲 ;wi 为 第 i 个 气象 因子 变量 ,本 文 包括 平均 风速 平均 气温 、 
平均 相对 湿度 和 日 照 时 数 。5, 的 正 负 反映 E7, 变 化 与 气象 因子 变化 的 一 致 性 ,5, 为 正 值 表 明 E7, 随 该 气象 因 
子 的 增加 而 增加 ;5, 绝 对 值 大 小 反映 E7, 的 变化 对 相应 气象 因子 变化 的 敏感 程度 ,绝对 值 大 表明 E7, 对 该 气象 
因子 变化 的 敏感 性 大 。 

气象 因子 对 E7, 变 化 的 贡献 量 计算 "1, 即将 单个 气象 因子 的 敏感 系数 与 该 因子 的 多 年 相对 变化 率 相 乘 ， 
贡献 量 最 大 的 气象 因子 为 引起 ET, 变 化 的 主导 因子 。 表 达 式 为 : 


EE™ 
日 Eee 
《 
[| 
A 
人 
OD 


S$, = lim 
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Con, = 5, RC, (4) 
Rd (5) 

式 中 ,Con, 为 气象 因子 vi 对 ET 变化 的 贡献 量 ;RC, 为 vi 的 多 年 相对 变化 率 ;av 为 60 年 vi 的 平均 值 ;Trend 为 
逐年 变化 率 ,由 趋势 分 析 法 计算 获得 。 
2.4 插值 方法 

本 研究 采用 ArcGIS 软件 中 协同 克 里 格 插值 方法 (Co-Kriging Interpolation ) ,地 理 数据 选用 空间 分 辩 率 为 
90 mx90 m 的 毛乌素 沙 地 数字 高 程 模型 (SRTM3 DEM ) 数据, 将 其 作为 地 形 因 子 参数 ( 协 变量 ) 对 研究 区 ET， 
及 其 对 气象 因子 的 敏感 系数 进行 空间 搬 值 。 搬 值 计 算 中 协 变异 函数 表达 式 为 : 


说 = 二 三 区 Ce) = Zs + [Zw) ~ Zr + 和] (6) 
式 中 ,yj() 为 协 变异 函数 ;NV'(h) 为 分 离 距 离 有 的 Z,(x) 和 Z(x) 的 样本 对 数 ; 
协同 克 里 格 法 的 插值 表达 式 为 : 
z" (wo) 三 Na) 二 Dal) (7) 


式 中 ,z*(%o) 为 xo 处 预测 变量 的 估算 值 ;z,(x,) 为 预测 变量 的 测量 值 ;z,(%) 为 协 变量 的 测量 值 ;A ,A 为 分 别 
为 z,/ 和 zz, 的 权重 ;m,g 为 分 别 为 预测 变量 和 协 变量 测量 点 数量 。 

对 空间 插值 结果 的 精度 和 误差 采用 t 检验 方法 进行 统计 分 析 ,P 值 大 于 0.05 ,未 达到 显著 差异 水 平 ,表明 
插值 估算 结果 与 实测 值 不 存在 显著 差异 ,协同 克 里 格 插值 结果 可 信 。 
2.5 数据 分 析 

本 文 利用 Matlab (R2012b) 统计 软件 ,应 用 Mann-Kendall 趋势 检验 法 对 E7T, 和 气象 因子 变化 趋势 进行 显 

著 性 检验 ;应 用 敏感 性 分 析 方 法 量化 分 析 主 要 气象 因子 (平均 风速 .平均 气温 、 相 对 湿度 和 日 照 时 数 ) 对 ET, 变 

化 的 影响 程度 及 其 贡献 量 ; 应 用 协同 克 里 格 插值 方法 对 研究 区 ET, 及 其 对 气象 因子 的 敏感 系数 进行 空间 插 
值 ,分 析 有 7 及 其 对 气象 因子 的 敏感 系数 空间 分 布 特征 。 


3 结果 与 分 析 


3.1 毛乌素 沙 地 参考 作物 节 散 量 及 气象 因子 时 间 变 化 趋势 
3.1.1 年 值 变化 趋势 

由 表 1 可知 ,毛乌素 沙 地 ET 多 年 均值 为 1048.81mm,ET, 最 大 值 在 定 边 站 (1242.35mm,2006 年 ) ,最 小 值 
在 横山 站 (851.64mm,1959 年 ) ,逐年 变化 率 为 0.54 mm。Mann-Kendall 趋势 检验 结果 表明 ,毛乌素 沙 地 87, 呈 
缓慢 上 升 趋势 ,其 中 , 定 边 和 鄂尔多斯 站 E7T, 上 升 趋势 达到 显著 性 水 平 ;而 鄂 托 克 底 和 伊 金 霍 洛 旗 站 多 年 FT， 
呈现 显著 下 降 趋势 。 

图 1 为 近 60 年 毛乌素 沙 地 年 87, 及 主要 气象 因子 的 变化 趋势 。 距 平 曲线 显示 多 年 57, 呈 波动 变化 过 程 ， 
ET, 在 1955 一 1968 年 呈现 下 降 趋势 ,1969 一 1975 年 稍 有 上 升 ,1976 一 1996 年 呈现 下 降 趋 势 ,1997 一 2006 年 逐 
渐 上 升 ,2007 年 以 后 又 呈现 下 降 趋势 。 近 60 年 内 毛乌素 沙 地 ET 的 波动 范围 在 906.82 一 1142.37mm 之 间 , 多 
年 相对 变化 率 为 3.11%, 随 时 间 呈 现 不 显著 增加 趋势 。 相 对 湿度 呈现 不 显著 变化 ,多 年 平均 值 为 50.91% ,多 
年 相对 变化 率 较 小 ,为 -0.71% 。 风 速 和 日 照 时 数 随时 间 整 体 呈 现 显著 降低 趋势 ,多 年 平均 值 分 别 为 2.72m/s 
和 7.94h, 距 平 曲线 显示 1983 年 以 后 风速 和 日 照 时 数 年 均值 明显 低 于 两 者 多 年 平均 值 ,多 年 相对 变化 率 分 别 
为 一 17.62% 和 一 6.05%。 人 气温 随时 间 呈 现 显 著 上 升 趋势 ,多 年 平均 值 为 7.75%C ,多 年 相对 变化 率 为 25.55%。 
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表 1 毛乌素 沙 地 各 站 点 1955 一 2014 年 参考 作物 年 蒸 散 量 特征 值 及 变化 趋势 


Table 1 Characteristic value and variation of annual reference crop evapotranspiration in Mu Us desert during 19SS 一 2014 


参考 作物 蒸 散 量 画 ee 人 ge . 
2 定 边 县 。 鄂尔多斯 市 。 哇 托 克 旗 。 横山 县 。 盐池 县 伊 金 震 洛 旗 。 榆林 市 。 毛乌素 沙 地 
Evapotranspiration 
平均 信 
1103.85 1043.44 1069.01 1043.89 1058.53 1050.84 1014.05 1048.81 

Mean value/ mm 
最 小 值 

971.76 893.80 945.88 851.64 915.31 952.50 891.85 906.82 
Minimum value/ mm 
l= 下 
最 大 . 1242.35 1187.72 1214.03 1214.86 1173.38 1167.67 1142.67 1142.37 
Maximum value/ mm 
ZE 责 
变化 率 3.44 0.81” =122.** 0.84 0.0048 一 1.87 0.84 0.54 
Rate of change/ mm 


* ,* 站 分 别 表示 变化 趋势 通过 了 置信 和 度 95% 和 99% 的 显著 性 检验 
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1 1955 一 2014 年 毛乌素 沙 地 参考 作物 蒸 散 量 及 气象 因子 年 值 变化 趋势 


Fig.l Change trend of annual reference crop evapotranspiration and meteorological factors in Mu Us desert during 19SS 一 2014 


3.1.2 ”季节 变化 趋势 

毛乌素 沙 地 E7T, 及 主要 气象 因子 具有 明显 的 季节 动态 变化 特征 ( 表 2)。E7T, 夏 季 最 高 ,多 年 平均 值 为 
5.23 mm/d, 冬 季 最 低 ,为 0.70 mm/d。 春 、 冬 两 季 多 年 平均 ET7, 星 现 显 车 上 升 趋势 ,逐年 变化 率 分 别 为 0.0036 
mm/d 和 0.0019 mmxd , 而 夏 、 秋 两 季 多 年 平均 ET, 变化 不 显著 。 风 速 春 季 最 高 ,为 3.34m/s, 冬 季 最 低 为 
2.47m/s。 相 对 湿度 秋季 最 高 为 57.08% ,春季 最 低 为 39.81%。 气 温和 日 照 时 数 季节 变化 同 ET, 一 致 。Mann- 
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Kendall 趋势 检验 结果 表明 ,四 季风 速 均 呈 现 显著 下 降 趋 势 ,气温 呈现 显著 上 升 趋势 。 春 . 冬 两 季 相 对 湿度 呈 
现 显著 下 降 趋 势 。 夏 、 秋 冬 三 季 日 照 时 数 显著 下 降 。 
3.2 ”毛乌素 沙 地 参考 作物 蒸 散 量 的 空间 分 布 

从 图 2 可 以 看 出 近 60 年 毛乌素 沙 地 平均 有 7 空间 分 布 基本 呈现 自 西向 东 递 减 , 地 区 性 差异 明显 ,西部 高 
蒸 散 量 区 ( 哪 托 克 旗 .盐池 县 和 定 边 县 ) BT, 变化 范围 为 1058.53 一 1103.85 mm ,东部 低 蒸 散 量 区 (鄂尔多斯 市 、 
伊 金 霍 洛 旗 .榆林 市 .横山 县 和 靖 边 县 ) 变化 范围 为 1014.05 一 1043.44 mm。 不 同 季节 ET, 的 空间 分 布 与 年 E7， 
空间 分 布 基本 一 致 ,呈现 西部 较 高 ,东部 较 低 。 

表 2 1955 一 2014 年 毛乌素 沙 地 参考 作物 蒸 散 量 及 气象 因子 季节 变化 趋势 


Table 2 Change trend of seasonal reference crop evapotranspiration and meteorological factors in Mu Us desert during 1955—2014 
春季 Spring 夏季 Summer 秋季 Autumn 冬季 Winter 


平均 信 ”变化 率 ”平均 信 变化 率 ”平均 值 。 ”变化 率 ”平均 值 变化 率 

参 藻 散 量 /( mm/d 

过 作 物 8 Ca) 3.62 0.0036“ 3.23 -0.0016 2.00 0.00047 0.70 0.0019™* 
Evapotranspiration 
平均 风速 /( m/s 

I Me 3.34 一 0.012 一 2.83 一 0.0058 一 2.48 一 0.0061 一 2.47 一 0.010 一 
Mean wind speed 
平均 气温 / 民 C 

均 飞 温 9.28 0.033 一 21.34 0.019 一 7.71 0.035 一 一 7.28 0.0553 一 
Mean temperature 
平均 相对 湿度 /4 

3 加 a 39.81 —0.086” 55.63 —0.038 57.08 —0.018 50.65 一 0.071” 
Mean relative humidity 
日 照 时 数 Solar duration/h 8.40 0.00039 8.81 一 0.013 7.47 一 0.011 6.92 -0.012™* 


* ,# #* 分 别 表示 变化 趋势 通过 了 置信 和 度 95% 和 99% 的 显著 性 检验 


3.3 ”参考 作物 蒸 散 量 与 气象 因子 的 敏感 性 分 析 

3.3.1 参考 作物 蒸 散 量 对 气象 因子 的 敏感 系数 A 
由 表 3 可 知 ,E7, 年 变化 对 风速 .气温 和 日 照 时 数 的 

敏感 系数 为 正 值 ,这 说 明 其 他 气象 因子 不 变 的 条 件 下 ， 3 

E7, 将 随 着 以 上 气象 因子 的 升 高 而 增加 ,而 E7v 对 相对 


湿度 的 敏感 系数 为 负 值 ,表明 相对 湿度 的 升 高 将 抑制 。 38 参考 作物 闲散 量 
E87, 增加 。 通 过 比较 敏感 系数 绝对 值 可 知 ,E7, 变 化 对 Es 
各 气象 因子 的 敏感 程度 依次 为 风速 > 日 照 时 数 > 气温 > ;。 i 
相对 湿度 。 EL 
不 同 季节 E7, 变 化 对 气象 因子 的 敏感 性 与 全 年 有 2 
所 差异 。 春 、 秋 两 季 E7, 变 化 对 日 照 时 数 最 为 敏感 ,其 2 毛乌素 沙 地 参考 作物 燕 散 量 空间 分 布 
次 为 相 对 湿度 夏 冬 两 季 ET 变 人 化 对 相 对 湿度 最 为 敏 Fig.2 Spatial distribution of reference crop evapotranspiration in 
感 ,其 次 为 平均 风速 。 Mu Us desert 
表 3 毛乌素 沙 地 不 同 季节 气象 因子 敏感 系数 
Table 3 Change trend of seasonal sensitivity coefficients for meteorological variable in Mu Us desert 

季节 平均 风速 平均 气温 平均 相对 湿度 日 照 时 数 

Seasons Mean wind speed Mean temperature Mean relative humidity Solar duration 

春季 0.25 0.27 -0.43 0.93 

夏季 0.77 0.62 -0.82 0.39 

秋季 0.37 0.37 =0.71 0.81 

冬季 0.79 0.63 一 0.95 0.44 

全 年 Annual 0.69 0.68 -0.62 0.68 


http://www.ecologica.cn 


201704.00154v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


1972 生 态 学 报 37 卷 


毛乌素 沙 地 E7, 对 气象 因子 的 敏感 性 存在 区 域 差 异 。 如 图 3 所 示 , 整 个 区 域内 ,7 对 风速、 气温 和 日 昭 
时 数 的 敏感 系数 均 为 正 值 ,而 对 相对 湿度 的 敏感 系数 为 负 值 。 空 间 分 布 上 ,风速 和 气温 敏感 系数 的 低 值 区 分 
布 在 毛乌素 沙 地 的 西北 部 地 区 ,同时 该 地 区 也 是 相对 湿度 和 日 照 时 数 敏感 系数 的 高 值 区 ,这 说 明 该 地 区 ET 
的 变化 对 风速 .气温 和 相对 湿度 敏感 程度 小 ,而 对 日 照 时 数 的 敏感 程度 大 。 风 速 敏感 系数 的 高 值 区 和 相对 湿 
度 的 低 值 区 同时 分 布 于 毛乌素 南部 地 区 。 和 气温 敏感 系数 高 值 区 位 于 毛 马 素 沙 地 的 东南 部 地 区 ,日 照 时 数 敏 感 
系数 的 低 值 区 分 布 于 毛乌素 沙 地 的 西南 部 。 


N 
平均 风速 
高 : 0.45 
低 : 0.01 
0 100 200km 
1 一 上 
| | | | | 
1l07"E 108°E 109%E 110E 111oE 107E 108°E 109%E 110E 111°E 
40°N 大 
39° 
平均 相对 湿度 时 
38° - 敏感 系数 敏感 系数 
高 : -0.31 高 : 1.70 
一 低 : -1.43 一 低 : 0.16 
37° 0 100 200km 200km 
| | | | | 
107。 108。 109。 110° 111°E 107° 108。 109° 110° 111°E 


3 ”毛乌素 沙 地 各 气象 因子 敏感 系数 空间 分 布 


Fig.3 Spatial distribution of sensitivity coefficients for each meteorological variable in Mu Us desert 


3.3.2 ”参考 作物 蒸 散 量变 化 的 主导 因子 

自然 条 件 下 ,气象 因子 对 ET, 的 影响 是 相互 联系 又 相互 制约 的 共同 决定 结果 ,气象 因子 对 有 7 变化 的 贡献 
量 能 够 定量 化 分 析 单个 气象 因子 对 ET 的 影响 ,57 变化 的 主导 因子 即 贡献 量 最 大 的 气象 因子 。 气 温 对 ET， 
年 变化 的 贡献 量 最 大 ,为 17.37% ,是 毛乌素 沙 地 PT, 变化 的 主导 因子 ,其 次 为 风速 (-12.15% ) .日 照 时 数 
(-5.11% ) 和 相对 湿度 (0.44% ) ,由 此 得 出 , 近 60 年 气温 的 显著 升 高 引起 £7, 的 增加 ,而 风速 和 日 照 时 数 的 显 
著 下 降 导 致 57, 减 少 。 

由 于 不 同 季 节 的 气候 特征 不 同 ,k7, 变 化 的 主导 因子 随 季 节 的 变化 发 生 改 变 。 夏 季风 速 对 E7, 的 贡献 量 
最 大 ,为 -8.34% ;春秋 、 冬 三 季 气 温 对 BT 的 贡献 量 最 大 ,分别 为 5.77% 、10.15% 和 27.70%。 由 此 可 知 , 春 、 冬 
两 季 PIT, 的 显著 增加 归 因 于 气温 的 显著 上 升 ;而 夏 、 秋 两 季 E7, 变 化 不 显著 ,是 由 于 风速 和 日 照 时 数 的 显著 下 
降 , 其 对 E7, 变 化 的 负责 献 量 与 气温 上 升 的 正 贡献 量 相抵 消 。 

空间 分 布 上 ,ET 变化 的 主导 因子 存在 区 域 差异 ,不 同 区 域内 气温 和 风速 对 7, 变化 的 贡献 量 均 较 大 ,其 
次 是 日 照 时 数 ,而 相对 湿度 的 贡献 量 最 小 。 气 温 作 为 主导 因子 影响 ET 变化 的 主要 区 域 集中 在 毛乌素 沙 地 东 
部 的 鄂尔多斯 市 .榆林 市 和 横山 县 ,以 及 西北 部 地 区 的 鄂 托 克 旗 , 风 速 作为 主导 因子 影响 £7, 变 化 的 主要 区 域 
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集中 在 毛乌素 沙 地 西南 部 的 定 边 县 和 盐池 县 ,以 及 东南 部 的 伊 金 霍 洛 旗 。 由 于 气象 因子 的 多 年 相对 变化 率 与 
敏感 系数 共同 决定 贡献 量 的 大 小 ,所 以 影响 ET 变化 主导 因子 的 空间 分 布 与 敏感 系数 分 布 有 所 差异 。 


4 讨论 


由 于 研究 区 域 及 时 间 段 不 同 ,西北 地 区 E7, 变 化 对 气象 因子 的 啊 应 程度 存在 区 域 差异 。 已 有 人 研究 表明 ， 
内 蒙古 额济纳 绿洲 £7, 变 化 对 太阳 辐射 敏感 程度 最 高 ” ,黄土 高 原 地 区 E7, 变 化 与 气温 和 日 照 时 数 极 显 著 相 
关 '” ,而 石 羊 河流 域 57, 变 化 与 相对 湿度 相关 性 最 好 "”" 。 本 文 基于 最 新 年 代 气象 数据 计算 敏感 系数 得 出 , 毛 
乌 素 沙 地 ET 变化 对 相对 湿度 敏感 程度 最 高 。 虽 然 敏感 系数 可 以 反映 E7, 对 气象 因子 变化 的 敏感 程度 ,但 是 
不 能 衡量 单个 气象 因子 对 ZE 变化 的 实际 贡献 水 平 , 即 定量 的 评价 各 气象 因子 对 E7, 的 影响 大 小 。 本 文通 过 
计算 气象 因子 对 87, 年 变化 的 贡献 量 , 得 到 气温 和 相对 湿度 对 毛乌素 沙 地 BT7, 变 化 具有 正 贡 献 作 用 ,风速 和 日 
照 时 数 对 ET 变化 具有 负 和 贡献 作用 ;气温 是 &7, 变 化 的 主导 因子 ,其 次 是 风速 ,这 与 曹 雯 等 ' 引 基于 1961 一 2009 
年 气象 资料 分 析 西 北 干旱 半 干 旱 区 结果 一 致 。1955 一 2014 年 毛乌素 沙 地 风速 和 日 照 时 数 呈 现 的 降低 趋势 抑 
制 了 ET 增加 ,但 是 对 ET, 变化 贡献 量 较 高 的 气温 呈现 显著 上 升 趋势 ,气温 对 FT, 增加 的 促进 作用 上 略 大 于 风速 
和 日 照 时 数 的 抑制 作用 ,相对 湿度 的 相对 变化 率 较 小 导致 其 对 7, 增加 的 贡献 率 不 大 ,因此 ,毛乌素 沙 地 FT 
整体 呈现 不 显著 上 升 趋势 。 

空间 上 ,由 于 气象 因子 对 ET 变化 的 贡献 量 及 气候 特征 存在 区 域 差 异 ,毛乌素 沙 地 ET 空间 分 布 不 均匀 ， 
总 体 呈 现 自 西向 东 递 减 趋势 ,加 之 毛乌素 沙 地 地 形 由 西北 向 东南 部 倾斜 ,这 将 加 重 毛乌素 沙 地 西北 部 地 区 的 
水 资源 短缺 压力 和 干旱 程度 。 

由 于 毛乌素 沙 地 区 域内 气象 站 点 有 限 ,本 文 基 于 可 获取 的 典型 气象 站 点 数据 代表 该 区 域 气象 特征 ,将 气 
象 因 子 作 为 主导 因子 分 析 E7, 的 时 空 变化 特征 ,但 在 干旱 半 干 旱地 区 ,由 风速 ,气温 和 相对 湿度 变异 引起 的 平 
流 、 大 气 环流 年 际 变化 下 垫 面 的 均一 性 以 及 人 类 活动 均 影 响 E7, 的 变化 ,可 是 这 些 因素 尚 不 能 定量 化 计算 分 
析 , 因 此 ,本 文 研 究 结果 存在 一 定 程度 的 局 限 性 。 参 考 作 物 蒸 散 量 作为 实际 蒸 散 量 的 理论 上 限 ,是 衡量 实际 鞘 
散 量 的 重要 参考 ,由 于 干旱 半 干 旱地 区 受到 严重 水 分 限制 ,开展 该 地 区 实际 蒸 散 量 和 参考 作物 燕 散 量 的 对 比 
分 析 ,能够 进一步 深入 和 完善 干旱 半 干 旱地 区 的 水 循环 研究 。 


S 结论 


(1) 近 60 年 毛乌素 沙 地 ET 整体 呈现 缓慢 上 升 趋势 ,多 年 平均 值 为 1048.81mm;E7T, 年 内 变化 夏季 最 高 ， 
多 年 平均 值 为 5.23 mm/d ,冬季 最 低 , 为 0.70 mm/d。 毛 乌 素 沙 地 E7v 空 间 分 布 整体 呈现 自 西向 东 递 减 趋势 。 

(2) 各 气象 因子 对 ET 的 敏感 性 分 析 中 ,ET, 对 风速 ,气温 和 日 照 时 数 为 正 敏感 ,对 相对 湿度 为 负 敏感 。 
ET 变化 对 各 气象 因子 的 敏感 程度 依次 为 风速 > 日 照 时 数 > 气温 > 相对 湿度 。 春 秋 两 季 E87, 变化 对 日 照 时 数 
敏感 程度 最 高 , 夏 . 冬 两 季 E7, 变 化 对 相对 湿度 敏感 程度 最 高 。 空 间 分 布 上 ,毛乌素 沙 地 东南 部 地 区 为 气温 敏 
感 系数 高 值 区 ,西北 部 地 区 为 相对 湿度 和 日 照 时 数 敏感 系数 高 值 区 ,南部 为 风速 敏感 系数 高 值 区 。 

(3) 气 温 是 影响 毛乌素 沙 地 E7, 年 变化 的 主导 因子 ,其 次 是 风速 日 照 时 数 和 相对 湿度 。 近 60 年 气温 的 
显著 升 高 引起 E77, 的 增加 ,而 风速 和 日 照 时 数 的 显著 下 降 导 致 E70, 减少。 夏季 ET, 变 化 的 主导 因子 是 风速 ; 
春秋 、 冬 三 季 主 导 因 子 是 气温 。 空 间 分 布 上 ,毛乌素 沙 地 西南 部 地 区 E7, 变 化 的 主导 因子 为 风速 ,东部 地 区 
主导 因子 为 气温 。 影 响 87, 变化 主导 因子 的 空间 分 布 与 敏感 系数 分 布 有 所 差异 。 
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